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Abstract. The stellar populations of galaxies contain a wealth of detailed information. From 
the youngest, most massive stars, to almost invisible remnants, the history of star formation 
is encoded in the stars that make up a galaxy. Extracting some, or all, of this information has 
long been a goal of stellar population studies. This was achieved in the last couple of decades 
and it is now a routine task, which forms a crucial ingredient in much of observational galaxy 
evolution, from our Galaxy out to the most distant systems found. In many of these domains we 
are now limited not by sample size, but by systematic uncertainties and this will increasingly 
be the case in the future. 

The aim of this review is to outline the challenges faced by stellar population studies in the 
coming decade within the context of upcoming observational facilities. I will highlight the need 
to better understand the near-IR spectral range and outline the difficulties presented by less 
well understood phases of stellar evolution such as thermally pulsing AGB stars, horizontal 
branch stars and the very first stars. The influence of rotation and binarity on stellar population 
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, modelling is also briefly discussed. 
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+3 ', The luminous output of normal galaxies is ultimately generated by stars in various 
stages of stellar evolution — in isolation a trite observation, but this simple fact gives us 
i— —i, the opportunity to extract a vast amount of information from observations of galaxies 
through the modelling of their stellar populations. 

While the study of stellar populations in galaxies started with Baade's (1944) idcnti- 
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1. Introduction 



fication of two populations of stars in M32 and NGC 205, a rigourous study of the topic 
only commenced in the late 60's and early 70's ( Tinslev 1968) : Faberlll972|: ISe 
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l/"") ■ Il973f ) with Tinsley's Fundamentals of Cosmic Physics article (|Tinslevlll980D particularly 
influential. Since that time, the number of articles discussing stellar populations has 
f^) ■ risen rapidly so that today about 12% of all articles in the major journals mention stellar 
populations in their abstracts (Figure 1). 

The majority of this growth has been made possible through the development of simple 
models for the evolution of stellar populations that have found widespread use in a wide 
range of astronomical studies, from stellar clusters in the Milky Way to the most distant 
galaxies in the Universe. I will repeatedly refer to stellar population models below, and 
5—1 ' use ^ ms term loosely to refer to any model that predict the observational properties of any 
ensemble of stars. These start from models of stellar evolution (e.g. Bertelli et al 1994; 
see contribution by Cassisi these proceedings). By applying empirical colour corrections 
(e.g. Lejeune et al 1997) they can be placed on an observational Hertzsprung-Russell 
diagram. For a given star formation history and initial mass function (IMF), this can be 
sampled, either to produce a Monte Carlo realisation of an observed colour-magnitude 
diagrams (see Tolstoy et al 2009), or by convolution to create integrated properties of a 
stellar population (e.g. Fioc & Rocca-Volmerange 1997; Leitherer et al 1999; Vazdekis 
1999; Bruzual & Chariot 2003; Maraston 2005; Kotulla et al 2009). Finally, these models 
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Figure 1. Left panel: The fraction of papers in A&A, AJ, ApJ and MNRAS each year that 
mention "Stellar populations" (solid line) or "Population synthesis" (dashed line, scaled by a 
factor of 8) in their abstracts. Currently about 12% of all papers mention "Stellar populations" 
and 1.5% "Population synthesis". Middle panel: The mass-metallicity relat ion for SPSS DR6 
(cf. Tremonti et al 2004 for DR2) with a metallicity adjustment following !Brinchman n et aD 
( 2008). Right panel: The relation between stellar mass and the 4000A break in the SDSS DR7. 
This is an updated version of the influential study by iKauffmann et all (|2003bf ) , showing the 
bi-modality of the galaxy population. 



are compared to observations using a number of techniques, in itself an interesting topic, 
but for this review the details of this step are not crucial. 

In this contribution I will begin by briefly reviewing some of the achievements in extra- 
galactic astronomy made possible by our understanding of stellar populations, before 
I turn to the challenges for stellar populations studies in the coming decade. Given 
the format this will not be an exhaustive review, see the other contributions in these 
proceedings for more in-depth discussion of many of the topics. 



2. A series of successes 



2.1. Resolved stellar populations — extracting star formation histories 

From the very beginning, the study of objects where individual stars can be resolved has 
been central to stellar population modelling. The reason for this is that one can construct 
a colour-magnitude diagram (CMD) directly from the observations and this is sufficiently 
close to the theoretical Hertzsprung-Russell diagram that insight into the properties of 
the population of stars can be had fairly straightforwardly. 

By comparing a theoretical CMD to the observational one, it is possible to infer a 
number of prop erties fo r the stellar system being studied. The first such study was 
arguably that bv lMaederl (|l974l ) but the full power of the approach has only been realised 
in the last couple of decades, with an early application of this technique to infer the star 
formation history of Local Group dwarfs bv lTosi et all (|l99lh . 

The field blossomed through the use of Hubble Space Telescope (HST) to obtain very 
deep CMDs of dwarf galaxies throughout the Local Group (see Tolstoy, Hill & Tosi 2009 
for a revie w) and through the d evelopment of sophisticated algorithms for the analysis of 
CMDs (see lGallart et al.ll2005l) . This wealth of data has provided us with an impressive 
insight into the history of star formation in small galaxies in the Local Group and it is 
now clear that they show an enormous range in star formation histories. 
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2.2. The estimation of stellar masses and the mass assembly of the Universe 

There is now a vast body of deep photometric observations of the skjfj] covering a large 
stretch of cosmic time. By matching these observations to stellar population models, 



it is possible to infer the stellar mass of the distant objects 


(e.g. iGiallongo et al. 


1998 


Brinchmann & Ellis 2000: Bell & de Jong 2001 


; Drorv et al. 


2001c iFontana et al. 


2003 


Kauffmann et al.ll2003atlDickinson et al.ll2003; Maraston et al 


2006). This has led to stel- 



lar mass becoming the most important independent variable for galaxy evolution studies 
Notable app lications of this is the qu antification of a transition mass in the local galaxy 
popu lation ( Kauffmann et alJ,l2003bl ) and the mass-metallicity relation ( Tremonti et al 



2004) , updated versions of the key plots from those studies is shown in Figure [TJ 



When the stellar masses of galaxies have been calculated for a statistically well-defined 
sample, it is possible to calculate the stellar mass density in galaxies. By combining 
samples over a wi de range in redshift, we can infer the stellar mass assembly history of 
the Universe (e.g. | Brinchmann &: EllisH2000t iDickinson et al.ll2003c IFontana et al.l 120031 : 



Arnouts et al. 20071 . see Wilkins et al 2008 for a compilation). 



This has provided a very visual image of the assembly of mass, and with the large 
samples currently available it is possible to do this in bins o f mass and hence deriv e the 
evolution of the stellar mass function with cosmic time (e.g. iMarchesini et aL 2009t ) and 
we now know that the massive end of the stellar mass function was in place at much 
higher redshift than was expected a decade ago. 



3. The future — context 

Major progress has come in extra-galactic research in the rece nt decade through l arge- 
scale surveys of the sky, such as the 2d F Galaxy Redshi ft Survey (IColless et al . 2001) and 



the Sloan Digital Sky Survey fS DSS. lYork et all 2000) in the nearby Unive rse, and the 
Vimo s Very Deep Survey (VVDS, iLe Fevre et al.ll2004f l and DEEP-2 survey (jPavis et al 



2003) out to z ~ 1. The coming years will see a continuation of this trend with, for 



insta nce, the commencement of the various VISTA and VST surveys (see lArnaboldi et~al 



20071) . and several Dark Energy surveys (e.g. BOSSj, Wiggle^ and HETDEX[||). These 



developments are very exciting and it is clear that stellar population modelling will be 
key to the exploitation of the data coming out of these studies. 

However, while these large-scale studies undoubtedly will be very valuable for a wide 
range of scientific questions, they do not present stellar populations with very different 
challenges from what we already face and hence I will not focus on these surveys here. 
The situation is somewhat different with the Large Synoptic Survey Telescop^HI (LSST), 
which, if built, will provide repeat imaging of the northern sky — opening up the possi- 
bility for a systematic use of variability in the study of stellar populations. This is clearly 
a challenge for the coming decade but not one I will discuss further here. Neither will 
I discuss the opportunities offered by new long-wavelength facilities such as Herschel, 
ALMA or SKA, or high energy facilities such as Fermi and IXO, even though these all 
offer novel challenges an d in the future we w ill hopefully see them integrated into stellar 
population analysis (c.f. da Cunha et al. 20081 ). 



Instead I will focus on the major change in studies of the extra-galactic Universe 
that will be ushered in with the launch of the James Webb Space Telescope (JWST, 

f For a compilation see http://www.strw.leidenuniv.nl/-jarle/DeepFields/ 
X http://cosmology.lbl.gov/BOSS/ 
http : //wigglez . swin.edu.au/Welcome .html 
|| http: / /www. as .utexas . edu/hetdex/ 
ft http://www.lsst.org/lsst 
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Gardner et al. 20061 ) and the building of on e or more Extremely Large Telescopes (ELTs) 



such as the European ELtQ, the TMlH and the GMT^. These facilities will present a 
significant shift in emphasis relative to much of the stellar population modelling carried 
out at present. 

The JWST will be a truly transformative influence on astronomy and will firmly shift 
the focus of extra-galactic observations to longer wavelengths than is commonly used 
today. At wavelengths A > 2/zm JWST will be two orders of magnitude more sensitive 
than current facilities, opening up an entirely new area of wavelength space for extra- 
galactic studies. 

The proposed ELTs will have two advantages relative to today's lOm-class telescopes: 
Much higher sensitivity and much higher spatial resolution. The former will enable us to 
obtain much larger samples of spectra of distant galaxies and detect much fainter and 
more distant galaxies than we currently can. The improved resolution will be achieved 
with the help of adaptive optics systems working in the near- and mid-IR, and this will 
enable resolved stellar population studies to reconstruct star formation histories for more 
distant galaxies and in more crowded regions. 

Thus, while extra-galactic research in the coming decade will make use of much larger 
samples and the routine observations of very faint objects, the main gain relative to today 
will come from studies in the near- and mid-IR; both because JWST will operate at those 
wavelengths, but also because those are the wavelengths where adaptive optics, and hence 
improved spatial resolution, will work best. The challenge for stellar population studies 
will be to exploit these data in an optimal manner. I should emphasise that this does 
not mean that optical astronomy will become obsolete, merely that I don't view that 
wavelength range as presenting as many new challenges as longer wavelengths and hence 
of less importance for this review, although it will be discussed below. 



4. Challenges — or areas of interest 

4.1. Error estimates on stellar population models 

Models for stellar populations have reached a high degree of sophistication and ease of 
use. As mentioned above, this has led many researchers to use them to help interpret 
their observational data. With this wide-spread adoption comes an increased need to 
quantify the uncertainties in the models as well — stellar population models are rather 
involved and it is difficult for a user to assess the reliability of a particular prediction 
without input from the model builders. These error sources are many: Uncertainties 
in the stellar evolution models for normal stars, such as the treatment of turbulent 
mixing; uncertainties in atomic data; treatment of uncertain or intractable phases of 
stellar evolution such as horizontal branch (HB) stars, thermally pulsing asymptotic giant 
branch (TP-AGB) stars, post-AGB stars or various phases of binary evolution; stellar 
wind loss; uncertainties in stellar atmosphere calculations; mismatches between stars in 
observational libraries and the theoretical tracks they are matched to and observational 
uncertainties in the empirical data included in the models — just to mention a few. 

These model uncertainties are likely to dominate the error budget in analyses of high 
S /N galaxy spectra, and until we understand them, we will not be able to take full advan- 
tage of the best data a telescope can deliver. Thus a major challenge for model builders 
in the coming decade is to construct models that incorporate uncertainty estimates, or 

\X http: //www. eso . org/sci/f acilities/eelt/ 
http://www.tmt.org 
|| http://www.gmto.org 
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probably more realistically, indications of reliability for different predi ctions. Notable 
first attempts at this has been recently published by ( Conrov et al. 20091 ). who explored 



the consequences of uncertainties in out understanding of HB and TP-AGB stars on the 
predictions of population synthesis models, and Percival & Salaris (2009), who explored 
the impact of uncertainties in the calibration of fundamental stellar parameters. These 
are promising first steps but are still far from a comprehensive study. This is likely to 
improve in the future but an accurate treatment of uncertainties is likely to remain a 
major challenge for years to come. 

4.2. The near-IR — can we properly make use of future observations? 

We emphasised above that adaptive optics on ELTs and in particular the flight of JWST 
will ensure that much of the focus for future astronomical observations will move into 
the near- and mid-IR, both for photometry and spectroscopy. It is therefore essential to 
ensure that models accurately reproduce the rest-frame near- to mid-IR fluxes of stellar 
populations, and to understand better what near-IR spectral features are of importance 
for studying stellar populations. 

We expect that resolved stellar population studies making use of the much improved 
resolution offered by adaptive optics systems on ELTs will be widespread in the future. 
It is therefore paramount that we start carrying out similar studies with current facilities 
to learn how to make optimal use of ELTs for this kind of work. The results from such 
efforts have only re cently started to app ear. An examp le of this is the study of Galactic 
globular clusters bv lOriglia et"aH (|2008l NGC6440) and lMoretti et alj (|2009l NGC 6388) 



using NACO on theVLT, and the studies by Fiorentino et al (in prep.) of NGC 1928 and 
ICampbell et al.1 (2008) of R136 using the Multi-Conjugate Adaptive Optics Demonstrator 



(MAD, iMarchetti et al.ll2007l ) on the VLT. While very useful, these latter two studies in 
the Magellanic clouds do not yet probe down to the main-sequence turn-off and hence 
do not put very strong constraints on models, hopefully this can be remedied in the near 
future with deeper observations. 

Resolved stellar population studies in nearby galaxies will not only be carried out using 
photometric observations, but also using adaptive optics allied to near-IR IFUs, such as 
the proposed EAGLE for the E-ELT, IRS/IRMOS on the TMT and GMTIFS on the 
GMT. To make optimal use of these facilities it is clearly necessary to make use of spectral 
features in the near-IR. This is an area that has seen si gnificant effort recently, with in - 



depth studies of the CO band-head at 2.3^m by e.g. Marmol-Queralto et al.l ( 20081 ). 



and assembly of ne ar-IR atlases of stellar spectra such as the IRTF Spectral Library 
(jRavner et al. 2009). At the moment these studies cover relatively restricted ranges in 



stellar parameter space and we need a considerable effort both on the theoretical and 
observational side in the coming decade to ensure that we can make optimal use of 
cutting-edge facilities on ELTs. 

Studies of resolved stellar populations can often ignore problematic stellar phases by 
excluding them from analysis. For the majority of galaxies, which are unresolved, there is 
no such option. Thus fitting unresolved stellar populations in the near-IR can be expected 
to be challengin g. This is bor n e out by experience, although for old stellar populations 
(ellipticals/SOs), [c^rterelalJ (|2009h showed that a wide range of population synthesis 



models give consistent results when applied to broad-band photometry inc luding near-IR. 
In gen eral, the situation is more complex as was shown, for instance, bv lEminian et aTl 
(2008) in their comparison between models and UKIDSS + SDSS photometry. They 



pointed out that while some combinations of near-IR and optical filters could be well-fit 
by models, others could not. 

To focus our discussion further, let us consider the determination of stellar masses. As 
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emphasised above, this is an essential ingredient of most galaxy evolution studies today 
and it is therefore important to understand what systematic uncertainties can influence 
their determinations. The estimation of stellar masses fundamentally involve a fit to 
either broad-band photometry or spectral features to estimate M^/Lx in some band, X, 
and finally this is scaled by Lx to give the stellar mass. Any systematic uncertainties in 
the models might lead to biases in the estimates of M*/Lx- 

The most obvious systematic uncertainty is probably the IMF, but we will ignore 
this here due t o spa ce limitations. Focusing instead on the population synthesis models, 
Conrov et all (|2009h showed that systematic uncertainties in stellar population models 
probably limit the accuracy of stellar mass determinations to about 0.3 dex. In addition 
to this we know that different algorithms for how star f ormation histories a re treated 
can lead to systematic differences of the order of 0.2 dex (jPozzetti et al.ll2007l ). but that 
when the algorithms are the same, different models give the same stellar mass to within 
0.15 dex (based on a comparison of the mass estimates from Bruzual & Chariot 2003 
models (BC03) and Maraston 2005 models (M05) by Tojeiro et al 2009) . 

Many studies looking at the impact of including rest-frame near-IR data in the fit have 
focused on higher redshi ft galaxies. They generally find that the situation is less good. 



van der Wei et al 



(2006) compared dynamical and stellar mass estimates at both low 
and high redshift and found that mo dels did not g i ve ac ceptable results when near-IR 
data was included. More recently, e.g. iMuzzin et al. ( 2009h found that model predictions 
in the near-IR were the limiting systematic uncertainty in their study of z ~ 2.3 galaxies. 
Some of this could be due to the way very young stars are treated in the models as Conroy 
et al (2009) showed that there is a glaring discrepancy between the observed V — K colour 
of LMC clusters and the predictions of models at young ages, but we note that this relies 
age measurements from a broad range of sources with rather different age estimation 
techniques and dust corrections are non-trivial. A careful re-examination of this issue 
might be worthwhile. However it is generally accepted that part of this discrepancy in 
stellar mass estimates is likely to be due to the treatment of TP-AGB stars in the models. 
These stars dominate the near-IR luminosity of stellar populations with an age of a few 
Gyr. Since this is the typical age of galaxies at z ~ 2-3, it follows that uncertainties in the 
treatment of this population in the models can have severe impact on the interpretation 
of observational data. This was emphasised by Maraston et al. ( 20061 ) who highlighted 
the differences in the treatment of TP-AGB stars in M05 and BC03. Subsequent work 
has done much to improve the situation and Marigo & Girardi have recently published 
impro ved evolutionary tracks including a si gnificantly improved treatment of TP-AGB 
stars ( Marigo et al.ll2008t iMarigo fc Girardilf2007i see also the contribution by Margio in 
these proceedings). Despite this, it is fair to say that there is still considerable uncertainty 
associated to the treatment of TP-AGB stars at low metallicity and high redshift. It 
is difficult to see that major progress can be made with resolved stellar populations 
alone, rather one might have to combine resolved stellar population studies with "almost 
resol ved" studies using surface bright ness fluctuation studies in the near-IR and optical 
(e.g. lRaimondol[2009l : IlLee et al.lliooa ). 



4.3. Non-solar abundance ratios — interpreting high S/N data 

While most population synthesis models assume scaled-solar abundance ratios for their 
stellar tracks and atmospheres, it has long been recognised that variations in the abun- 
dance ratios of elements away from solar ratios can strongly affect any inferences one 
makes based on spectral features. Deviations from scaled-solar abundance ratios reflect 
a variation in past star formati on histories, in particular the relative contributions by 
Supernova Type II and la (e.g. Worthev et al. 1992 ; Trager et al. 20001 : Thomas et al 
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2003). This i s normally seen where the t ypical time-scale for star f ormation is shorter th an 
~ lGyr (e.g. Matteucci fc Recchi 200lh . as in elliptical galaxies ( Thomas et al. 20051) or 
at z > 4. 

Treating non-solar abundance variations has long been a challenge for stellar popula- 
tion models, but the last decade, in particular, has seen great progress. There is now a 
fairly well-developed machin ery for dealing with a- variation within the Lick index system 
using fitting functions (e.g. Tripicco fc Bell 19951 : Trager et al. 20001 ; Thomas et al. 2003 : 
Tantalo et al.ll2007t ISchiavonll2007l) and this has been extensively used to study the prop- 



erties of massive galaxies i n the nearby Universe (e.g. iThomas et al.ll2005t I Graves et al 
2009allbl ; ISnTith et al.l[2009h . 



There has also been gre at progress in the calculation of theoretical isochrones for 
non-solar abundance ratios ( Pictrinfcrni et al. 2004; ICassisi et al. 2004; iPietrinferni et al 



20061: IWeiss et al.ll2007t iDotter et all 20071 and th ere has also been significant progress i 



the calculation of stellar atm osphere models fe.g. lCoelho et alj|2005t iMunari et al.ll2005f ). 



see 



Martins fc Coelhol (j2007l ) for an in-depth comparison of some of these libraries. Taken 



together this has allowed th e construction of libraries suited to po pulation synthesis of 
integrated populations (e.g. ICoelho et al. 2007 ; Percival et al. 20091 ). These models are 
now sufficie ntly mature that the y can begin to be used for the interpretation of galaxy 



spectra (e.g. Watcher et alj|200"9T) . 



Despite the great progress made, there are still some major issues outstanding. The 
first is that considerable work must be done to understand what detailed features in 
the theoretical libraries can be fully trusted — this needs careful testing against high- 
resolution spectroscopy. Related to this, and also of great importance and interest, is the 
question of whether it is sufficient to just vary the a-elements as a block, or whether 
the individual elements must be varied one by one. It has been argued that the latter 
is indeed necessary to optimally extract information from high-S/N spectra, and some 
progress has been made recently ( Potter et al.l 12007 ; Lee et al. 2009). This does expand 
the available parameter space significantly and should be applied in a careful manner. 



4.4. The first stars and galaxies 

One of the major challenges for stellar populations in general will be to interpret ob- 
servations of the very high-redshift Universe from JWST. It is clear that it will become 
necessary to understand what kind of objects we are likely to see, will they be truly 
zero metallicity, or will they have very small metallicities? And even within the zero 
metallicity class there is now tho ught to be two distinc t classes of stars, referred to as 
Pop III.l and Pop III. 2 stars by Tan fc McKed ( 2008 ). dependin g on whether or not 
their formation was a ffected by radiation by an earlier population (jMcKee fc Tanl 12008 ; 
Greif fc Bromm 20061 ). This in turn has important implications for the IMF of the first 



stars, and therefore on the strength of the Hen 1640A line which has often been suggested 
as a probe of Population III stars ( Tumlinson et al] 12001 , 2003: Sc haerer 2003 ). 

In addition to these uncertainties, it has also recently been realised that the halos 
where the first stars form have a very high dark matter density. If the dark matter is 
its own anti-particle, self-annihilation is a possible en ergy source for the first stars, a 
topic that has seen considerable interest recen tly (e.g. Spolvar et al. 20081 Freese et al 



20081 : Iocco et al. 2008 : Ripamonti et al.l [2009). It is not clear whether the effect of dark 
matter annihilation on the first stars will turn out to be crucial, or indeed observable. 
Howev er it has been argued that similar eff ects might be observable in the Galactic 
centre ([Scott et alj|2009t iFairbairn et al.ll20081 ). thus this is an area with potential future 
implications. 
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4.5. The rest-frame ultraviolet — understanding massive stars 

In any galaxy with on-going star formatioan, the rest-frame ultraviolet (rest-UV) spec- 
trum is dominated by O and B stars. This can be viewed as both a disadvantage — our 
understanding of massive stars, whether their evolution is strongly influenced by rotation 
(e.g. Meynet et al 2009), or significantly influenced by binary evolution (e.g. Eldridge et al 
2008), is certainly lacking in some aspects. On the other hand it is an advantage because 
it offers us a direct probe of massive stars, providing us with a wealth of information to 
test and develop models. With the upcoming demise of HST, the coming decade will see 
a dearth of UV-sensitive facilities, but optical spectrographs on ELTs will provide very 
large samples of optical spectra of galaxies at z ~ 2-3, sampling the rest-UV. 

Thus it is a major challenge for stellar population models to be able to extract in- 
formation from these spectra (e.g. Rix et al 2004; Maraston et al 2009). As discussed 
by Leitherer (these proceedings), stellar atmosphere models for hot, massive stars are 
now reaching maturity and will likely be sufficient for most future studies. On the other 
hand, the evolution of massive stars, their mass-loss prescriptions and disentangling the 
influence of rotation and of binaries remain a challenge. 

It is interesting in this context to note that the light-gathering power of a 42m-class tele- 
scope is about 25 times larger than that of an 8m-class telescope (for resolved sources). 
Thi s is similar to the amplification seen in gravitationally lens e d Lyman-break galax- 



ies (Allam et al 



Lin et al. 2009: 



■2007t IBelokurov et all 120071 : ISmail et alj|2007t iBelokurov et alj l2009t 
Kubo et al.ll2009h . Thus by observing these objects, it is possible to 



test models on very similar data to that which w ill be obtained rou t inely with an ELT , 
something that seve ral groups have realised (e.g. iPettini et al ] l2000t iQuider et al.ll2009l : 
Hainline et al. 2009f ). This should clearly continue and a concerted effort for a survey of 



the nearby Universe with COS would also be very useful. 



4.6. The importance of binaries 
Most stars with mass > IMq are born in binaries and the binary fraction appears to 



increase with stellar mass (e.g. lLadall2006l ). thus there is no doubt that binary evolution 
must be important for understanding the stellar populations of galaxies at some level (see 
the contribution by Vanbeveren in these proceedings). However the inclusion of binaries 
into population synthesis models is problematic as it leads to a much increased parameter 
space and complexity. 

Including binary evolution is however crucial to understand some rare, but bright, 
evolutionary phases. This includes certain variable star phases, X-ray binaries and stellar 
remnants, most of which should ideally be included in population synthesis modelling. 
However in the UV-optical region, the main concern is with horizontal branch stars (see 
Moni Bidin et al 2008 for a review). These stars dominate the near-UV light in old stellar 
populations and understanding whether they are due to binary evolution (e.g. Han et al 
2007)or not has important implications for the analysis of elliptical galaxies. 

Binary evolution might also be crucial for the development of Wolf-Rayet stars, par- 
ticularly at low metallicities. Brinchmann et al (2008) showed that large number of 
Wolf-Rayet stars at low metallicity is a significant challenge to stellar evolution models. 
Eldridge & Stanway (2009) recently showed that this tension can be resolved by binary 
evolution effects, although similar results can be achieved by inclusion of rotation in the 
evolution of massive stars (e.g. Meyner & Maeder 2005). The relative importance of these 
two effects is not yet understood and is an important challenge for the future. 
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5. Conclusions 

Rather than a traditional conclusion, I would here like to end by summarising the 
challenges above into a set of rest wavelength ranges, highlighting the issues and some 
suggestions as to what might be useful studies to do — note that most of these suggestions 
are merely reinforcing already existing studies! 

Far-UV — Massive stars — what are realistic rotation and binary parameters. To 
do: Observations of lensed Lyman-break galaxies, COS observations of nearby stars and 
star-forming regions, more in-depth studies of stellar mass-loss 

Near-UV — Binary evolution, mass loss on the RGB. To do: COS and GALEX spec- 
troscopy of UV-upturn galaxies, inclusion of horizontal branch uncertainties in models. 
Combination of optical, near-UV and X-ray data in analysis. 

Optical — High-precision predictions. Non-solar abundance variations and impact 
on spectra, emission lines. To do: Careful and extensive testing of theoretical spectra, 
comparison to detailed spectroscopic analysis of nearby systems: what abundances can 
be reliably extracted from medium resolution spectra? 

Near-IR — Evolution of and spectra of luminous, cold stars. TP-AGB stars, RGs and 
RSGs, spectral features in the near-IR. To do: Further observations of resolved stellar 
populations in the near-IR with and without adaptive optics with careful comparison 
to results including optical data. Calibration spectral features in the near-IR. Obtain 
ACS optical data to back up future studies in near-IR with ELTs. In depth analysis 
of uncertainties in population synthesis models and careful comparisons to data in the 
nearby and distant Universe. 
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